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In existing literatures about the top-hat spherical collapse model of galaxy clusters formation 
in cosmology containing dark energies, dark energies are usually assumed not to cluster on this 
scale. But all these literatures ignored the current describing the flowing of dark energies outside 
the clusters which should exist under this assumption, so the conclusions of these literatures are 
worth further explorations. In this paper we study this model in QCDM or Phantom-CDM cos- 
mologies(flat) by assuming that dark energies will cluster synchronously with ordinary matters on 
the scale of galaxy clusters so the dark energy current flowing outside the clusters does not exist 
at all and find that in this case, the key parameters of the model exhibit rather non-trivial and 
remarkable dependence on the equation of state coefficients of dark energies. We then apply the re- 
sults in Press-Scheter theory and calculate the number density of galaxy clusters and its evolutions. 
We find that this two quantities are both affected exponentially by the equation of state coefficients 
of dark energies. We leave the study of this model with the assumption that dark energies do not 
cluster on the scale of galaxy clusters at all as the topic of another paper where similar conclusions 
will be obtained also. 



I. INTRODUCTION 

When Press-Schechter theory is used to study the num- 
ber density of galaxy clusters and its evolutions, there are 
two key parameters 0, needi ng t o be extracted from 
the spherical collapse models. 0,^3 are two early works 
studying this model. They considered it in a background 
universe without dark energies. [Tl EHl I2I I2I 
are some recent works containing discussions about this 
model in ACDM or QCDM cosmologies. In these exist- 
ing works, when dark energy is included, it is assumed 
that it does not cluster on the scale of galaxy clusters. 

However, if dark energy does not cluster on this scale, 
then when we write down the Einstein equations to de- 
scribe the evolution of an over-dense region which will 
develop into the future galaxy clusters, we should include 
a dark energy current in the energy momentum tensor to 
describe the dark energy component flowing outside the 
over-dense region. But, in all these existing literatures, 
such a current was ignored. To see the effects of such a 
current, there are two methods. The first is, adding it in 
the energy momentum tensor and solving the resulting 
Einstein equations. The second is, assuming that such 
a current does not exist at all and solving the appropri- 
ate Einstein equations to see the effects of dark energy's 
clustering on the scale of galaxy clusters. 

The purpose of this paper is to study this problem 
by the second method. We leave the study of the first 
way as the topic of another paper |28|. in that paper we 
will get similar results as in this paper. After a sim- 
ple introduction to Press-Scheter theory in section^ we 
study the spherical collapse model in section IIIII in de- 
tails. Our study shows that when dark energy is assumed 
to cluster synchronously with the ordinary matters so the 
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dark energy current does not exist on the scale of galaxy 
clusters, the key parameters of the model exhibit rather 
non-trivial and remarkable dependence on the equation 
of state coefficients of dark energies. We then apply the 
results obtained in this section in Press-Scheter theory 
and calculate the number density of galaxy clusters and 
its evolutions in section IIVI We find that, under this 
assumptions, the number density of galaxy clusters as 
well as its evolutions are both affected exponentially by 
the equation of state coefficients of dark energies. We 
then study the possibility of determining the equation 
of coefficients of dark energies by measuring the number 
density of galaxy clusters. We end the paper in section 
[V] which contains the main conclusions of this paper and 
some discussions. 

Our study is performed for time-independent equation 
of state coefficients of dark energies w which ranges in 
—00 ~ 0. But in some cases we will use terminology 
Quintessence or Phantoms regardless how much w is. 

II. PRESS-SCHECHTER THEORY 

Press-Schechter Theory ^ij predicts that the fraction 
of volume which has collapsed at a certain red-shift z is 
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f co u(M(R),z) = 2 / dSe- s2 ^ 2 ^\ 

y / 2na(R, z) Js c 
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Here, R is the radius over which the density field has been 
smoothed, cr(R, z) is the rms of the smoothed density 
field Q . 8 C is the threshold of density contrast at time z 
(count stopping time) beyond which objects collapse. 
By the notations of 0, chapter 7, 

cr 2 (i?,z) = 

2 [Jk^T 2 (k)W 2 (k-R)dk] + 
CT s k n »+ 2 T 2 (k)W 2 [k-8h- 1 Mpc)dk] £>f(l) U 



where n s is the spectral index of perturbation powers, 
n a = 1 corresponds to scale invariant power; T(k) is the 
transfer function, in this paper we will use the BBKS fit- 
ting formulaes for it is the smooth- 
ing window function, we will use the top hat window in 
this paper; + is the growth function of linear 
perturbation theory, 



£>i(«) 



H{a) I da'[a'H{a')}- 3 . (3) 
o 



Operationally 21], S c is defined as the extrapolation of 
primordial perturbations to the collapse epoch using the 
growth law of linear perturbation theory, i. e., 
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Pmc(a) 

Prnb(a) ~' -Dl(o) 



Di(a c ), 
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where p mc and p m b are the matter densities of clus- 
ter and background respectively. It can be shown that 
Di(a) a ^o — ► a. Using the method of we can show 
that 



= (1 — a ■ a), 



(5) 



where a p is the scale factor of the cluster, we use the 
normalization for a p so that when a — > 0, a p ~ a, see 
section Ull Al Substituting cq© into eqJ3J we get 



3a • Dx(a c ). 



(6) 



By studying spherical collapse model we will give an an- 
alytical formulae for the calculation of a. 

We comment here that, because Di(a) a ^o - * a, as 
long as the definition eq|0J is to be of any sense, we must 
have [y^ — l]a->0 ~ a - So eq(0 holds regardless what- 
ever the background cosmology model is. Physically, this 
is because, any realistic cosmological model behaves like 
a flat matter dominated cosmology at early times, so they 
have the same limit behavior as a — > 0. Mathematically 
eq(|5j is only derived out explicitly in the OCDM cos- 
mologies 0| and ACDM cosmologies [25| and j2||. 

It should be noted besides the partition of the cos- 
mological components and the observational epoch, the 
value of S c is also quite dependent on the choice of 
smoothing window used to obtain the dispersion a(R, z) 
We will not consider this effect in this paper. 



a constant curvature parameter to describe the evolution 
of the over-dense region. However, if this is the case, 
then when we write down the Einstein equation = 
SnGT^, we should include a dark energy current -Q 
current- in the energy momentum tensor to describe the 
flowing of dark energy outside the over-dense region. But 
and HI did not do so. 
To see the effects of such a current, we have two 
choices. The first is directly including such a current in 
the energy momentum tensor and solving the resulting 
Einstein equations, we will do this way in [2^|. The sec- 
ond is assuming that such a Q-current does not exist and 
solving the resulting Einstein equations to see the effects 
of Quintessence's clustering on the properties of galaxy 
clusters. In this method, Quintessence will cluster or 
expand synchronously with matters in the over-dense re- 
gion, just as we usually do in the case of universes. This 
is an indirect but easier method. We will take this way 
in this paper. 

Considering an uniform over-dense region in the back 
ground of a flat cosmology containing general Quin- 
tessence component, the background satisfies 
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For the over-dense region, starting from the metric 
ansaltz 

ds 2 = -dt 2 + U{t,r)dr 2 +V{t,r)(dd 2 +sm 2 9d<p 2 ), (8) 
using Einstein equations, we can prove that 
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where A; is a constant independent of time while a p sat- 
isfy: 
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a p d p + 2d 2 + 2k = 4TrG[p mc + p Qc - PQ c }a p . 

The difference of cas ljlOj l and yields 

r d Dl9 k 8irG , . 
[-? + -r = -5-PW- 



(10) 
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III. SPHERICAL COLLAPSE MODEL IN QCDM 
COSMOLOGIES 



A. Einstein Equations and Energy Conservation 

In [23| and j^, it is stated that the energy hold by 
the dark energy-Quintessence inside the the over-dense 
region varies independently of the region's radius, and 
the curvature parameter of the over-dense region is time 
dependent, so we have no Friedman- like equation, with 



In eq(7|), a denotes the scale factor of the background 
universe, p m b and pQb denote the density of matter and 
Quintessence in it respectively. While in eal|5|)- (|12|) . a p 
denotes the scale factor of the over-dense region, the 
future galaxy clusters. p mc , pq c denote the density of 
matter and quintessence inside the over-dense region re- 
spectively, while pQ c denotes the pressure produced by 
the Quintessence. By the statement of |23| and j23|, the 
curvature parameter k of the over-dense region is time- 
dependent. This is equivalent to say, the metric func- 
tion U(t,r) cannot be factorized as eq@. If this is the 
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case, then not only Freidmann equation (I12|l . but also 
Raychaudhuri-equation H10JI . or eq(A2) of 23] will follow 
the usual way. 

In practice, we need taking a normalization for the 
definition of scale factor a p of the over-dense region. We 
choose the convention so that a p ~ a when a — *■ 0. Un- 
der this convention, the value of k in eq© is not 1 defini- 
tively, its value depends on the total energy density inside 
the over-dense region. But, it is a constant independent 
of time. 

Using the metric ansaltz eqlJEJ and 10 , the energy con- 
servation equation Tff = can be integrated to give 



Pn 



a constant, 



(13) 



where w is the equation of state coefficients of the 
Quintessence(— 1 < w < 0) or phantoms(u; < —1). 
In this paper, we only consider Quintessence/phantoms 
with time-independent w models. During the formation 
of galaxy clusters, the mass of the system is conserved 
at least approximately. This directly leads to the con- 
clusion that, the density of Quintessence inside the over- 
dense region does not vary independently of the radius 
of the region. If the equation of state coefficients sat- 
isfy — 1 < w < 0, then as a p decreases, pq c increases, 
Quintessence clusters. If w < —1, then as a p decreases, 
Pq c also decreases, i.e., phantom is repulsed outside the 
over-dense region, it anti-clusters. Contrary to our rea- 
sonings here, in |27j . the energy held by Quintessence in 
the over-dense region does not conserve, see eq(14) of it. 

In [29l I30L l3lj . it is pointed out that on cosmolog- 
ical scales, a smoothly distributed, time-varying com- 
ponent violates equivalence principle, so is un-physical. 
Here, on the scale of galaxy clusters, we see that to pre- 
serve Einstein equations, we should either include a Q- 
current in the energy momentum tensor, or assume that 
Quintessence (or phantom, even time independent vac- 
cum energy) clusters (anti-clusters) just like matters. As 
is well known, when the total universe is looked as a 
cluster of size Hubble scale, such a current which de- 
scribes Q-component flowing outside the Hubble horizon 
is usually assumed not exist. Galaxy clusters are a class 
of objects whose formation and evolution is still in the 
linear perturbation theory applicable area and in the en- 
ergy composition of these objects, dark energy occupies 
a rather big part, just like our biggest cluster - our to- 
tal universe. So the assumption that Quintessence will 
cluster on the scale of galaxy clusters like ordinary mat- 
ters is at least partly the case. We will also see in j2^| 
that, if dark energy is assumed not to cluster on the this 
scale, then all the effects displayed in this paper will be 
enforced instead of weakened. 



region will expand as the background universe expands; 
but because this region's over-dense, its expanding rate 
will decrease and stop doing so at some middle times; 
then it starts to shrink because of self-gravitating, the 
final fate of this over-dense region is a singular point. 
But in practice, when this region shrinks to some degree, 
the pressures originate from the random moving of par- 
ticles inside the over-dense region will balance the self- 
gravitation and the system enters the virialization period. 
In theoretical studies, it is usually assumed that the viri- 
alization point is coincident with the collapse point of the 
ideal model on the time axis. 

According to Press-Shceter theory eqQ), if an over- 
dense region is to be virializated at some time a c , its 
density-contrast should be no less than 8 c {w,£l m o,a c ). 
While to relate the mass of a galaxy cluster with 
its characteristic X-ray temperature, the ratio of clus- 
ter/background matter densities at the virialization point 
is a very important parameter, 



A c (w,n m0 ,a c ) 



Pmc,c 
Prnb.c 



(14) 



To see the effects of Quintessence/Phantom's 
clustering/anti-clustering on the formation and proper- 
ties of the galaxy clusters, we would like to calculate this 
two important parameters in this section. To calculate 
S c (w, !] ra oi dc), starting from eqQ and lfT2T) . let 
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Pmc,ta 
Pmb,ta 
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x = — ,y = 



PQb,ta 



Pmb,ta 

and using the fact that [a p ]ta = we can get 

dy 3 = Cy-' + K 1+w y- x - 3w - (C + »C 1+W ) 



(15) 
(16) 

(17) 



Then using eq(j3J) and the approximate mass conserving 
condition we can prove 



and 
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^- [j/s— o,x#o] = C 3 (1 ~ 2a • a ta x) 



(18) 



(19) 



Substituting eqs©, (JTSJl and JTJJJ into eq{T7|), taking 
limit x — > and preserving only the first order small 
quantities in x, we can get 



B. The Key Parameters of the Spherical Collapse 
Model of Galaxy Clusters Formation 

In the ideal model, if there is an over-dense region in 
a flat background universe, then at very early times, this 



Substituting this result into eq© we can get 
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S c (w,n m0 ,a c ) = gO^tC* + ^ ■ C i+w ]Di(a c ) (21) 
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Now let us come to A^,s calculation. According to the 
definition ea (|14|l . 

a i f \ \ Pmc,c Pmctadp ta a p,c / - a p,ta & c 

A c (w, n m0 , a c ) = - — - = p ; _ 3 = C 



Pmb,i 



3 -3 
Pmb,taQ>ta@'C 



a p,c a ?a 



(22) 



To calculate the second factor of the above equations' 
right-most part, we can use energy conserving condition 
and virial theorem. Assuming that at the collapse point, 
the system virialized fully, then we have the following 
relations: 



inetic.c 



■yUmm.c + UmQ.c— -Uqq. c (23) 
1, 



+ 2C/ m g iC + -Uqq. c 

= U mm j a + U m Q t t a + Uqqjo, (24) 

where U mm , U m Q and Uqq denote the matter-matter, 
matter-Quintessence and Quintessence-Quintessence 
gravitation potentials respectively. The subscripts ;C 
and ,4a mean that the quantities carrying them take 
values at the collapse and turn around time respectively. 
Explicitly, ea l|24|) can be written out as: 



mc,c + ^Pmc,cPQc,c + PQ c ,c) a p,c 
2 i „ „ i 2 



^(Pmc^ta 

+ Pmc,taPQc,ta + PQc,ta) a p,ta ( 25 ) 

dividing each side of this equation by the relation 
J ta and making a little deformation 
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we get 



^mc,ta p 
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Pmc,c 



Prac.c 



= 2(1 



PQc,ta 
Pmc : ta 



PQc.ta s 

n 2 

rmc.ta 



(26) 

Using energy conservation law and the approximate mass 
conserving condition, we can change ea (|26[l into the fol- 
lowing form: 



|- a P.*a jQ | ^j- a y.t« j3w | t2r a P,ta-\6w\ _ ■>, , , , , ,2, 



j ) = 2(i+e+a 



(27) 



where 



£ = vC (28) 

Numerical studying shows that the third term on both 
sides of ea(|2T)l can be neglected, just as [23 did. But, for 
a general valued w, even neglecting this two terms, this 
equation cannot be solved analytically. So we preserve 
this two terms in the equation and resort to numerical 
methods to calculate ^£^2-. 

a P ,a 

From eqs(j2U and i|22|) . (|27[) we see that, given w, £l m o 
and a c , to calculate S c and A c , we need to express ata, v 
and C as functions of w, f2 m o and a c . For v 

1 Qmb,ta 



a, 



;b,ta 



1 



1 + ^ 



1 — "mO 



(29) 
(30) 



While for at a , we can use the fact that t c — 2t ta and 
cqQ to set up an integration equation 



da 



1 + VQd 



r-3w 



-O-ta 

= 2 1 da' 



1 + v a'- 3w 



(31) 



In the case of w = — 1, ea ()31|l can be solved analytically: 

,-,1/3 

(32) 



ata 



+ v {) a\ - 1 
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While in the case of general valued w, ea l|31|l should be 
solved numerically to get the function at a (a c ). 

For C's calculation, there is two method. The simplest 
is directly integrate ea l|17|l to set up an equation: 



dx\ 



1 + vx 



dy[ 



-3w -I 



(33) 



Solving this equation will tell us £'s dependence on w 
and v. Combining the results with ea l|29|) will gives us 
C's dependence on w, U m o and a c . The second method 
is to use cqs{7Jl, ijTUt and the notations in ea (|T5)l - ljPo)| to 
set up an eigen- value problem 



[*]-~o = i- ■ C* 

x 

y\ x= i = 1, y'\ x =i = 0. 
d 2 y _ dy 1 A + dfl mb 1 
efe 2 dx 2 x dx f2 



(34) 



1 



+ o[(i + 3#] 



^|-3(l+ll))^l+w ^(j: ^mb) _|_ C-E^mfcj _ q 



(35) 



Solving this eigen-value problem can also tell us £'s de- 
pendence on w and v |4|. For a given value of to, we 
can get £'s dependence on fl rn b,ta by solving ea (|33H or 
l|34 |l -l(3*5 )l and fitting the results as ( — ((Q m b,ta)- For 
example, in the case of w = — 1, 



a _ r37r l2o -0.7384+0.245in„ 
^ — i~\ ^ l mb,ta 



(36) 



But for the general cases when £l m b,ta and w both vary, 
we cannot fit the results into a simple function, see Fig^ 
It should be noted that this kind of fitting formula ca l|36|) 
only coincide with numerical result when O m ^ ta is not 
too small, for example, Sl m b,ta > 0.05. 



C. Numerical Results, Effects of w on the Key 
Parameters of the Model 



We provide our numerical results in FIG I1I3I From 
FIG^we see that, when the dark energy's clustering is 
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FIG. 1: f 's dependence on Q m b,ta and ui. Left panel, consider- 
ing Quintessence's clustering effects. It is solved from eg (13 3H 
or l|35jl . Right panel, the results of assuming Quintessence 
does not cluster on the scale of galaxy clusters but neglect the 
Q-c: urrent flowing outside the clusters, solved from eq(A9) of 



A c (w,n m0 ,a c =l) A c (w,n m0 ,a c =l) 




FIG. 2: A c 's dependence on Slmo, w and a c . The upper right 
part is solved from eq(5)-(6) of |23| where Quintessence is 
assumed not to cluster on the scale of galaxy clusters but 
neglect the Q-current flowing outside the over-dense region. 
The others are solved from eg 1331 or 1351 of this paper. 

considered, £ depends on the value of w remarkably. For 
a fixed £l m b,ta, C nrs t increases then decreases as w de- 
creases in the range [—1.7, -0.41, this forms a strong con- 
trast with the result of [23 , which is also depicted in the 
figure. In FIG0 besides similar dependence of A c on w 
and £! m o, we can still see that for a given f2 m o and w, 
A c decreases as a c decreases, i.e. the early an over-dense 
region virialized, the less is its matter density over that 
of the background. This is easy to understand, because 
the earlier an over-dense region virializes, the less dark 
energy is contained in the background cosmology, so the 
weaker the anti-cluster effects will be and the less dense 
a region is required to keep balance by self-gravitations. 

In FIG|21 we display <5 c 's dependence on w, £l m o and 
a c . Unlike in FIG0 we did not compare S c 's depen- 
dence on w, f2 m o in this figure between the cases where 
Quintessence is assumed to cluster and the contrary. By 
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FIG. 3: S c 's dependence on w, Q m o and a c . Quintessence's 
clustering effects are considered. 

the method of [23|, eq(A12)-(A13), it is difficult to get the 
precise value of 5 c (w, fi m o, a c ). Since this method needs 
to take the value of [^] at some point x = xq — > 0, xq 7^ 0, 
while the coefficients of the differential equation which 
is used to solve £ has singularities near this region, the 
quantity S c (x) calculated by this method depends on the 
selection of xq rather artificially. By the way, we note 
that eq(A12) may missed a term £ on the right hand 
side. By the statement of [23j, 8 C ~ 1.6 — 1.686, where 
— 1 < w < 0. Here, in our FIG|21 we see that, if a c = 1, S c 
may be as large as 2.1 when w « —1.7. As a c decreases, 
S c also decreases and goes to the limit value 1.686, what 
ever value w and Q m o takes. This is trivial, because the 
earlier an over-dense region virializes, the more likely is 
the back ground universe a totaly matter dominated one. 

IV. THE NUMBER DENSITY OF GALAXY 
CLUSTERS AND ITS EVOLUTIONS 

A. Theoretical Formulaes 

According to Press-Schechter theory, the comoving 
number density of clusters which have collapsed (i.e., viri- 
alized) at certain red-shift z and have masses in the range 
M - M + dM could be calculated: 

„ ( M, 2 ),M = -^i|A^ p[ _^ ldM 37, 

in this equation: ptot is the energy density of background 
universe, M — ^-R 3 ptot, so the factor denotes the 
average number density of clusters with mass M. The 
other factors in this equations are just obtained by dif- 
ferentiating eqQ with respect to M. In the original pa- 
per of Press and Schechter 14], the quantity p and M 
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appearing in ea (|37|l only refer to density and mass of 
matters. But under our assumptions, Quintessence and 
matters cluster synchronously so we have to understand 
them as the total density and mass. 

To relate the mass of a cluster with its characteris- 
tic X-ray temperature, consider a virialized over-dense 
spherical region in the background universe containing 
quintessences, according to virial theorem, we have 



kinetic, vir 
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(38) 
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where V 2 ir is the mean square velocity of particles in the 
cluster when the system is fully virialized and a p is the 
physical radius of the cluster. From this equation, we 
have 
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y vir ^ 
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FIG. 4: Effects of w on the number density v.s. temperature 
function of galaxy clusters when z — 0. The larger is as or 
the larger the function value will be. All four figures have 
S2 m0 = 0.27, h = 0.71 and n s = 1.0. 



Using relation: 



2 

vir 



(39) 



where k B is the Boltzmann constant, m p is the mass of 
proton, while p,m p is the average mass of particles, [3 is 
the ratio of kinetic energy to temperature. So we have 
mass-temperature relation: 



M = 
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GH{z) 



5[3k B T 1 
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H 2 



[im p f(z) 
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H(z) 
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with 
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T2 Qb , c A c " J )A c ]3 x 



rnb.c + ^Qb,c^c 

H(z) = H [Sl m0 (l + zf + (1 - r» m0 )(l + z)^ 1+w ^ 

(42) 

and A c given by ea l|22|) . 

Just as p3l | pointed out, since the mass-temperature 
relation is red-shift dependent, simply substituting 
ea H41|) into ea l|37|) cannot give us correct number density 
of clusters in a given temperature range today. Instead, 



we should first find out the virialization rate and mul- 
tiply it by the mass-temperature relation then integrate 
over red-shift 



n(T, z)dT 
1 



Ptot d\na dlna S Cf S 2 . 
MTdlnR dz j6XP[ 



2a 2 



]dzdT 



(43) 



From eqs(g3l, (glj, J2U and © we can see that besides 
the constant |j, n(T,z) will also depend on the cosmo- 
logical parameters w, f2 m o, h, n s and the normalization 
<78 of the cosmic density fluctuations. In principle, if 
we can measure the number density v.s. temperature 
relation precisely enough, by numerical fittings, we can 
determine all these parameters or some of their special 
combinations simultaneously from observations. 



B. Numerical Results, Effects of w on the Number 
Density of Galaxy Clusters and Its Evolutions 

In FIGQ] we display the effects of w on the galaxy 
clusters number density v.s. temperature relation when 
z = for two values of 4. From this figure we can see 
that the value of w affects the number density of galaxies 
exponentially. Just as we explained in section llll Al since 
the author of [23| and (2^ did not include the Q-current 
in their energy moment tensor at the right hand side of 
Einstein equations, which should appear under their as- 
sumptions, the conclusions there about the effects of w 
on the number density v.s. temperature relation is worth 
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FIG. 5: Effects of w on the evolution of number density of 
massive galaxy clusters. The larger is og, the more remarkable 
the effects of w on the evolution will be. Both figures have 
n m0 = 0.27, h = 0.71 and n s = 1.0. 



further explorations. Although we do not directly correct 
their analysis by including such a current, we only assume 
that Quintessence will cluster on the scale of galaxy clus- 
ters so that such a current do not exist at all, our results 
here may from the contrary tell us the importance of such 
a current on the scale of galaxy clusters. 

Comparing the number density v.s. temperature rela- 
tions between the two values of ^ case, we can see that 
this parameter also affects the relation strongly. This 
parameter ^ has physical meaning of efficiency of energy 
transformation from thermal-dynamic to x-ray. So the 
lower is this efficiency, the less possibly a galaxy cluster 
will be observed if its mass is not big enough. While 
the larger is the clusters, the less the number density 
of them will be. In [2l| and [23 , the authors fit the 
number density v.s. temperature relations to numeri- 
cal simulations to normalize this parameter and get that 
^ pa 1.3 ~ 1.4. By the notations of these works, the sym- 
bol fp is used to denote the same relevant quantit y a s 
our ^, our f u = 1, see the notations under eq(4) of |23| . 
Considering the time dep endence of mass-temperature 
relation explained by |23j, the value of this parameter 
is about 0.944 according to j^. Since we assume that 
Quintessence clusters synchronously with ordinary mat- 
ters, while the Q-component could not make contribu- 
tions to the x-ray emission processes. We should expect 
an even lower value of 41 than 0.944. In the numerical 



analysis in the next subsection, we will set ^ = 0.75, but 
will point out qualitatively the way a larger or smaller 
value of ^ affects our conclusions. 

In FIG[S]we displayed the effects of w on the evolution 
of number density of massive galaxy clusters for one value 
of M . Note that 



nM>M 



n{M',z)dM' 



(44) 



From the figure we can also see that the clustering of 
Quintessence affects the evolution of the number density 
of massive galaxy clusters exponentially, and that the 
larger is as, the more remarkable the effects will be. 
In this figure, what should interests us mostly may be, 



if w is too greater than —1 , the number density of galaxy 
clusters would increase as we look back to the past, while 
if w is too less than —1, the contrary is the trend. And, 
the larger is as, the strong this trend will be. Physi- 
cally, this trend can be easily to understand. Because, 
the more w — > 0, the more weakly Quintessence anti- 
clusters, so the more usually emerging of galaxy clusters 
could take place which will lead to the decreasing of num- 
ber of galaxy clusters. On the contrary, if w < — 1, then 
the more w — ► — oo, the more strongly Quintessence anti- 
clusters, so the more often splitting of massive galaxy 
clusters will take place which will lead to the increasing 
of the number of galaxy clusters. According to the cur- 
rent observations as we look back to past, we see 
smaller and smaller number density of galaxy clusters. 
This may implies that the equation of state coefficients 
of dark energies should not be greater than — 1 too much. 

Note that, we do not directly consider the emerging 
and splitting of galaxy clusters. However, our explana- 
tions about the two contrary evolution trends of the num- 
ber density of massive clusters above is still valid. Be- 
cause, there are such configurations where initially "two" 
clusters are closely located, if w < —1, then as time goes 
on, this two clusters go far and far and finally really be- 
come two clusters. However, if —1 < w, then as time 
goes on, these "two" clusters will not separate until to- 
day and have to be counted as one. This difference of 
"2" — >2 and "2"— >1 processes can be sensed by Press- 
Scheter theory itself through the effects of w on the key 
parameter S c displayed in FIG|21 with no use of direct 
consideration of gala xy clusters emerging and splitting 

processes us nana. 



C. Constraints From Observational Results 

Although the clustering of dark energy on the scale of 
galaxy clusters is only a logically possible assumption. 
Since this assumption predicts that the number density 
of galaxy clusters and its evolutions are affected expo- 
nentially by the equation of state coefficient w of dark 
energy. If theoretical predictions are to be coincide with 
observations, we should have constraints on w, £l m o and 
as which may be rather different from that comes from 
other observation such as CMB ^(j or large scale matter 
power spectrum |ll| . Does these observation leave spaces 
for our assumption? 

Motivated by this fact, we try our best to fit the obser- 
vation results of Q (with a factor of 2 corrected by 0) 
and into theoretical formulaes to find constrains on 
w, fl m o and &s- Our fitting method is minimizing x 2 , for 
details please refer to p|. However, constrained by our 
computation powers, we do not vary all the parameters 
involved in this problem and look for the best compo- 
sition. Instead, we choose to fit the measurements into 
theoretical formulaes in three cases. In the first case, we 
focus our attentions on as and Q m o but fix the other pa- 
rameters. In the second case we focus our attention on as 
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and w, while in the third case we focus our attentions on 
flmO and w. We also tried best to fit the four parameters 
4, w, fl m o and as simultaneously, but the convergence 
of the fitting result is so poor that we can not get any 
meaningful conclusions so we will not provide numerical 
results in this case. 

In the data source we used in this paper, ||| is the 
observed number density v.s. temperature relation at 
z « 0.05. The data from this source only gives errors in 
the ^(number density v.s. temperature function values) 
axis. We will fit it with our theoretical formula ea (|43[) . Q 
is the observed number density v.s. red-shift relation of 
galaxy clusters whose comoving 1.5/i _1 Mpc radius inside 
mass is greater than 8 x 1O 14 /i _1 M0. The data from 
this source have large errors in both the x(red-shift) and 
y (number density) axis. We will neglect its errors on the 
x axis although it is very large. We also neglecte the fact 
that the number density provided by this source is only 
that of the clusters whose comoving 1.5/i _1 Mpc radius 
inside mass is greater than 8 x 10 14 /i _1 M Q and think 
that is the number density of clusters whose total mass 
is greater than 8 x 10 14 /i~ 1 Af Q . This approximation will 
make our fitted as greater than the real value. For other 
observational results in this area, please refer to and 
the references therein. 

For reasons explained in the above two paragraphs, we 
should not look at the numerical results of this subsection 
too seriously. What should be emphasized is the con- 
sistency between our assumption that dark energy will 
cluster synchronously with ordinary matters and the ob- 
servational results of WMAP and SDSS. 

We display our confidence level analysis of the best fit- 
ting in the first case in FIG|S] where except as and fi m Q, 
the other parameters are fixed as w = — 1, h = 0.71, 
n s = 1.0 and ^ = 0.75. In this figure there are two 
points which is worth noticing or could be criticized by 
peoples. The first is, even use only the observed number 
density v.s. temperature relation of Q, we determined 
the cosmological parameters as and O m o to some degree. 
This forms a strong contrast with the conclusion of [2l| 
and |23| . which says that only using this observational 
data we cannot determine any one of this two parame- 
ters but only a special composition of them. We find if 
we neglect the normalization factor of the growth factor 
-P i (1) in e all2*l) . we will qualitatively reproduce the results 
of [21| and [23, please see our FIG[7| We also find that to 
get the results of |2l|, Figure 1, we need to artificially set 
the initial value of either of the two parameters 6o or k, 
eq(A16) and (A21) of |2l|. The second important point 
which could be criticized by peoples is, our results indi- 
cate that as ~ 0.5, which is exceptional low comparing 
with the results of WMAPpH and SDSS[H. We think 
this is because when WMAP or SDSS make their best 
fittings, they did not consider the perturbation of dark 
energies. If the perturbation of dark energies are consid- 
ered, we expect WMAP and SDSS would also predict a 
lower as than their current reported value. We note that 
the current version (V4.5.1) of the program CMBFast 
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FIG. 6: ACDM cosmology best fitting of galaxy clusters' ob- 
servation. The black boxed region denote the 68.3% confi- 
dence region of best fittingthe number density v.s. tempera- 
ture relation reported by [y], the +ed region denote the 99% 
confidence region of the same fitting. The blank boxed denote 
the 68.3% confidence region of best fitting the number density 
v.s. red-shift relation reported by while xed region, the 
99% confidence. If ^ takes values greater than 0.75, the black 
boxed and -l-ed region moves down left. 



have been able to calculate the Cosmological Microwave 
Background Anisotropy (CMBA) spectrums when dark 
energies are perturbed , but WMAP and SDSS did not 
use this version of CMBFast to fit their results. 

As is well known, an exceptional low as implies a 
strongly biased universe, while as = 1 implies a non bi- 
ased universe, which has been observed by the galaxy 
astronomers. So our result indicate that, on the scale 
of galaxy clusters, the mass distribution is rather dif- 
ferent from that on the scale of galaxies. This is just 
what we should expect, because we assume that on the 
scale of galaxy clusters, Quintessence clusters like ordi- 
nary matters, while on the scale of galaxies, the clustering 
effects of Quintessence are usually neglected by peoples. 
Physically this is because on the scale of galaxy clus- 
ters, Quintessence contribute a rather large part to the 
total energy of our studying objects, while on the scale 
of galaxies, the contributions of Quintessence to the to- 
tal energy of the studying objects is really negligible. It 
should be noted that our conclusion as ~ 0.5 is obtained 
by fitting the observational result of the number density 
v.s. red-shift relation into theoretical formulaes and has 
no dependence on our choice of ^ = 0.5. 

In FIG 03 we display the confidence level analysis of our 
best fitting in the second case where except as and w all 
the other parameters are fixed as £l m o — 0.27, h — 0.71, 
n s = 1.0 and ^ = 0.75. About this figure we would like 





give two comments. The first is, our results give an ex- 
ceptional low value of as, the physical reason of which 
we have explained in the previous paragraph. The sec- 
ond point is, although w affects the number density v.s. 
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FIG. 7: Still ACDM cosmology best fitting of galaxy clus- 
ters' observation, but ignored the normalization denominator 
Di(l) of the growth factor appearing in eq@. The meaning 
of different symbols appearing in the figure is the same as that 
of FIG© 
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FIG. 8: QCDM cosmology best fitting of the galaxy clusters' 
observation, except w and aa, all the other parameters are 
fixed as fi m0 = 0.27, 4 = 0.75, h = 0.71 and n s = 1.0. The 
+ed region denote the 68.3% confidence region of best fitting 
the number density v.s. temperature relation reported by [6j, 
the blank boxed region denote the 68.3% confidence region of 
best fitting the number density v.s. red-shift relation reported 
byi. 



temperature functions and the number density v.s. red- 
shift relations of massive galaxy clusters exponentially, 
see FIG 0] and El requiring theoretical predictions to co- 
incide with observations does not give constraining on w. 
This is because, the effects of increasing w and decreas- 
ing as on the number density of galaxy clusters counted 
each other so this two parameters degenerate. 

In FIG03 we display our confidence level analysis of the 
best fitting in the third case where except w and f2 m o, 



FIG. 9: QCDM cosmology best fitting of the galaxy clusters' 
observation, except w and Q m o, all the other parameters are 
fixed as as = 0.5, | = 0.75, h = 0.71 and n s = 1.0. The 
black boxed region denote the 68.3% confidence region of best 
fitting the number density v.s. temperature relation reported 
by 13, the +ed region denote the 99% confidence region of 
the same fitting. The blank boxed region denote the 68.3% 
confidence region of best fitting the number density v.s. red- 
shift relation reported by [9J, while the xed region, the 99% 
confidence. If ^ takes values greater than 0.75, the black 
boxed and +ed region moves down left. 



all the other parameters are fixed as as = 0.5, h = 0.71, 
n s = 1.0 and % = 0.75. On the basis of these priors, 
we get w = -1.08 ± 0.09, Q m0 = 0.27 ± 0.09 to 99% 
confidence level. For reasons explained in the first two 
paragraph of this subsection, we should not look at this 
numerical results too seriously. 



V. CONCLUSIONS 

We studied the top-hat spherical collapse model 
of galaxy clusters formation in the flat QCDM or 
Phantom-CDM cosmologies under the assumption that 
Quintessence or Phantom clusters or anti-clusters like 
ordinary matters. We found that under this assump- 
tion, the key parameters of the model exhibit rather non- 
trivial and remarkable dependence on the equation of 
state coefficients w of Quintessence or Phantoms. We 
then applied the results in Press-Scheter theory and 
calculated the number density v.s. temperature func- 
tion and the evolution of the number density of massive 
galaxy clusters and found that these two Quantities are 
both affected by w exponentially 

For the number density v.s. temperature function of 
galaxy clusters, we found that the nearer w — > 0, the 
larger the function value will be, the more w — -» — oo, the 
smaller the function value. While for the evolution of 
the number density of massive galaxy clusters, we found 
that if w is too less than —1, the number density de- 
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creases as we look back to the past. On the contrary, 
if w is too greater than — 1, the number density of mas- 
sive galaxy clusters increases as we look back. Using this 
fact, we studied the possibility of determining w by the 
observational results of galaxy clusters. Constrained by 
our computation powers and the big errors of the obser- 
vational data, we cannot determine all the parameters 
involved in this problem independently. However, on the 
basis some priors, we get that to 99% confidence level, 
W = —1.08 ± 0.09. On the other hand, if we fix w = — 1 
as priors, we find that cr 8 s» 0.5 which is exceptional low 
comparing with the reported value of WMAP and SDSS. 
This is because WMAP and SDSS both did not consider 
the perturbation of dark energies. 

Since from the beginning, we assume that dark en- 
ergy will cluster synchronously with ordinary matters, we 
have no dark energy current flowing outside over-dense 
region galaxy clusters. Although this is too strong an as- 
sumption, this assumption simplifies our discussions and 
makes it a self-consistent one. As a comparison, we point 
out that the discussions in the previous literatures which 
assume that Quintessence does not cluster on the scale 
of galaxy clusters but neglect the Q-current flowing out- 
side the over-dense region is not self-consistent and the 
conclusions of these literatures are worth further explo- 
rations. Our studies do not correct this problem directly, 
but may from the contrary indicate the importance of 



such a current's effects on the number density of galaxy 
clusters and its evolutions. 

In [2^, we will dispose the assumption that dark en- 
ergy cluster synchronously with ordinary matters and 
will find that, in that case, the equation of state coeffi- 
cients of dark energy affects the number-density of galaxy 
cluster even more strongly! 

As a discussion, we would like to point out that the 
actual case should be, dark energy should cluster on the 
scale of galaxy clusters but not do so synchronously with 
ordinary matters. So that is a case lies between the syn- 
chronous clustering with ordinary matter and not clus- 
tering on the scale of galaxy clusters at all. What ever 
the actual way is, the equation of state coefficients should 
affect the number density of galaxy clusters remarkably. 
So measuring the number density of galaxy clusters and 
its evolutions is a potential effective way to determine 
the equation of state coefficients. 
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